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Abstract 

We show explicitly that supersymmetric E$ Grand Unified Theory with a Higgs 
sector consisting of {27 + 27 + 351/ + 351' + 78} fields provides a realistic scenario for 
symmetry breaking and fermion mass generation. While gauge symmetry breaking 
can be achieved without the 78 field, its presence is critical for a successful doublet- 
triplet mass splitting. The Yukawa sector of the model consists of only two symmetric 
matrices describing all of quark, lepton and neutrino masses and mixings. The fermion 
mass matrices are computed at low energy and a fit to the second and third generation 
masses and mixings is performed. We find a good numerical fit to the low-energy 
data. Thus, this model, having 11 superpotential parameters, alongside the two 
symmetric Yukawa matrices, seems to be the best realistic candidate for a minimal 
renormalizable supersymmetric Eq unified theory. 
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1 Introduction 

There are three types of unified groups (simple groups containing the Standard Model gauge 
group) that admit chiral representations for fermions: SU(1V) for N > 3 [Tj, SO(4A + 2) 
for N > 2 and E 6 [2J. The first two have been thoroughly studied in many versions, 
renormalizable as well as non-renormalizable, with and without supersymmetry. Strange 
enough, Eg, the only exceptional group which contains chiral representations, has been 
mainly ignored over the last forty years. Apart from a few exceptions where symmetry 
breaking was analyzed (via renormalizable potential with 78 and two 27s in Ref. [3] and 
via renormalizable superpotential with 78, 27 and 27 in Ref. 0), only the Yukawa sectors 
have been touched upon without using explicit constraints arising from symmetry breaking 

0 El El El El HD!. 

Recently, two of us have tried to fill this gap by suggesting a possible minimal renormal¬ 
izable supersymmetric Eg model mi. The Higgs sector is composed of pairs of fundamental 
27 + 27 and the two-index symmetric 351' + 351' representations. Much to our surprise we 
found that, although successful in the symmetry breaking pattern to the Standard Model 
(SM) gauge group, this sector is unable to provide two light Higgs doublets of the minimal 
supersymmetric Standard Model (MSSM), or technically said, to perform the doublet- 
triplet (DT) splitting. The impossibility of DT splitting was due to the identical nature 
of the expressions for a massless doublet and for a massless triplet; this result was ob¬ 
tained by an explicit computation, and we are not aware of a shortcut explanation of this 
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unusual feature]]] To put this result in perspective, consider the minimal renormalizable 
supersymmetric breaking sectors (able to break to the SM group) for other groups. In the 
minimal renormalizable SU(5) superpotential (consisting of the adjoint 24), there are no 
SM doublets. In the corresponding minimal SO(IO) case (made out of 210, 126 and 126) 
one has in principle enough doublets, but only one (in 126) is coupled to fermions. In a 
minimal Eg case with 27 + 27 + 351' + 351', there are many more doublets which couple to 
fermions, yet the naive DT splitting fails (a completely different approach to DT splitting 
in E 6 can be found in [12]). We thus proposed in HU to add another 27 + 27 pair. In 
this way, the Yukawa sector consisted of three Yukawa matrices: the matter fields (3 copies 
of 27 p) have one coupling to the 351' and a coupling to each of the two Higgs-like 27s. 
Though we did not perform it explicitly, we strongly believe that in this model a realistic 
fit of the Yukawa sector is possible to obtain, but the model is essentially not predictive 
due to the large number of free parameters in the three Yukawa matrices. 

The purpose of this paper is to find an E 6 theory with only two Yukawa matrices and 
thus a simpler Yukawa sector than the model in HU, while keeping the symmetry breaking 
sector as simple as possible. The idea is to add to the minimal breaking sector (instead of 
the extra 27 + 27 pair) another multiplet, which preferably has the following properties: 

1. It does not couple to the matter bilinears (and so is not a 27, 351 or 351'). 

2. It contributes to the symmetry breaking vacuum. 

3. It increases the size of the mass matrices for weak doublets and/or color triplets. 

The minimal such multiplet is the adjoint representation 78, for which all of the above 
properties hold true. We will see that its inclusion is enough to allow the doublet-triplet 
splitting in the theory, thus leading to a theory with the correct (SM) vacuum and the 
correct (MSSM) low energy field theory. 

The Yukawa sector in this model consists of two matrices only. Assuming 3 generations, 
the down quark and charged lepton sectors have 3 extra vector-like fields, while the neutrino 
sector has 3 extra vector-like lepton doublet-antidoublet pairs, as well as 6 SM singlets, 
all coming from the extra fields in the decomposition of the fundamental 27 into the SM 
subgroup. The projection of all these fields into the usual 3 light generations will be 
performed explicitly and shown, in the simplified case of only two generations, to provide 
a realistic fit of the masses and mixing angles. 

Special attention needs to be paid to the overall neutrino mass scale. In minimal 
SO(IO) [13j m EE] this scale poses a serious problem and eventually rules out the low- 
energy supersymmetry scenario [EH El HEl [EH]: a too low SU(2 ) r breaking is disallowed 
by unification constraints and/or D = 6 proton decay [20], while a too large breaking scale 
makes the right-handed neutrinos too heavy and thus the seesaw mechanism ineffective 
to produce a large enough scale for the light neutrinos. In our model this conclusion is 
avoidable because there are many possible fields that can break SU(2)^. While the neutrino 
mass scale can be adjusted by choosing the vacuum expectation value (VEV) of the (1, 3,10) 
of the Pati-Salam (PS) SU(2) L x SU(2) R x SU(4)^7, located in the 126 of SO(IO), which in 
turn is found in the 351' of Eg, all other fields’ masses depend also on several other VEVs. 
We can thus arrange a lighter VEV of this (1, 3,10) still avoiding dramatic changes in the 
renormalization group (RG) gauge running. 

lr The sparsity of SU(5) breaking vacuum expectation values in the respective mass matrices is likely to 
be the fundamental origin of this problem. 
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We arrange the paper in the following way: we define the model and specify the terms 
in the superpotential in section [2} We find a suitable vacuum solution, which breaks Eg into 
the SM group in section [3j perform the doublet-triplet splitting in section [4] and compute 
the mass matrices in the Yukawa sector and identify the low energy content of the theory 
in section [5| We then use the obtained results to perform a 2 generation fit of the masses 
and mixing angles in section [6] We finish with a discussion in section [7} Five different 
appendices collect various definitions and technical results. All explicit computation in Eg 
has been performed using methods from [2TJ and HU; also, [22j 23- [21] might also prove 
useful to the reader. 

For ease of use we stick to the following color convention in this paper: red denotes 
VEVs coming from the spontaneous symmetry breaking of the unified theory at the scale 
Mqut , while blue denotes VEVs coming from electroweak (EW) symmetry breaking at the 
scale rriEw- 


2 Defining the model 

The renormalizable Eq SUSY GUT that we consider here is an alternative to the model 
in EH- It is motivated by the fact that the minimal Higgs sector in a renormalizable SUSY 
Eg model, which can break Eg to the Standard Model, is 351'+ 351'+ 27 +27. This minimal 
breaking sector cannot accommodate doublet-triplet splitting, however, and therefore needs 
to be extended to get a realistic model. One possible extension is an addition of a 27 + 27 
pair, which was analyzed in EH, while an alternative, taken in this paper, is to extend it 
by the representation 78 instead. The Yukawa sector in the present model will consist of 
only two matrices, unlike that in EH, which has three matrices. 

Our renormalizable Eg model thus contains the following: 

• The “fermionic sector” of three copies of a chiral supermultiplet 27, denoted by 27^, 
with i — 1, 2, 3. Also, we assume a Z 2 matter parity, under which the 27 p are —1, and 
the remaining chiral superfields are +1. With this symmetry, the ansatz (27 p) = 0 
is consistent with the equations of motion, which we shall adopt. 

• The “breaking sector” consisting of 351' + 351' + 27 + 27 + 78. 

The model under consideration is supersymmetric. The problem of SUSY breaking will 
not be considered, since it is (usually) an orthogonal problem to breaking the gauge group. 
One can imagine however, that whatever the detailed mechanism of SUSY breaking might 
be, we get soft SUSY breaking terms at energies not far above tuewj such that we get 
automatic unification of gauge couplings at Mqut assuming no extra particle states up to 
the unification scale. The effective low energy theory of our model will thus be MSSM. The 
soft SUSY breaking terms do not give contributions to the fermion masses at tree level. A 
fit to fermion masses and mixings can then be done without specifying the details of SUSY 
breaking, although the RG evolution of these parameters does depend on the SUSY scale 
(which we assume is of order few TeV). 
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The full superpotential of our model is 


Wfull = 777,27 27 27 + 777351 / 351/ 351' + 777 78 78 78 

+ Ai 351' 3 + A 2 35F 3 + A 3 27 2 35I 7 + A 4 2 7 2 351' 

+ A 5 27 3 + A 6 27 3 + A 7 27 78 27 + A 8 35I 7 78 351' 

+ \Y£ 2T F 27 j F 27 + \YiL_ 27 l F 27 j F 35I 7 . (1) 

Note that the definitions of 351' and 351' used here are switched compared to Slansky’s 
defintions [2SJ. 


3 Spontaneous symmetry breaking 

3.1 Equations of motion 

The Higgs sector 351' + 351' + 27 + 27 + 78 forms a realistic Higgs sector, which is able 
to break the gauge group from E e to SU(3)c x SU(2)x x U(l)y. We provide just such a 
vacuum below. 

First, note the decompositions of 27, 78 and 351' under SO(IO) x U(l): 

27 = 16(1) + 10(—2) + 1(4), (2) 

78 = 45(0)+ 16(-3)+16(3)+ 1(0), (3) 

351' = 1(8) + 10(2) + 16(5) + 54(—4) + 126(2) + l44(-l). (4) 

The representation 351' contains 5 SM singlets, 3 of which are SU(5) singlets (in 1, 
16 and 126 of SO(IO)), and 2 are part of a 24 under SU(5) (in 54 and 144 of SO(IO)). 
Similarly, the representation 78 also has 5 singlets, with one being a 24 under SU(5) (the 
one in 45 of SO(IO)), while the remaining ones are singlets under SU(5) (the 1, 16 and 16 
of SO(IO), as well as another one in the 45). The Higgs sector 351' + 351' + 27 + 27 + 78 
therefore contains 5 + 5 + 2 + 2 + 5 = 19 singlets in total. We list their VEV^jJin Table [lj 

With the ansatz (27 F ) = 0, the breaking part of the superpotential is 

H = 77735!' ^351'x351' + m 2 7 127x27 + m 78 hsx7& 

+ Ai / 351 '3 + A 2 + A 3 I 27 2 X 35i7 + A 4 I-2f x351 , 

+ A 5 I27 3 + A 6 1^3 + A 7 / 27x78x27 + As 7 3 5i/ x78x 35l7- (5) 

Explicit computation yields the following expressions for the superpotential invariants 

1 Notice that the Standard Model singlet VEVs have been denoted by u\, m, v, w and y. The notation 

from [4j is changed due to the states now being those which have well defined transformation properties 

under the SU(5) and SO(IO) subgroups of E 6 . The connection between the two notations is U\ oc a 1; 

v>2 oc tt9, w oc U3 — 63, v oc — -\/3a3 + 2a 4 — V 3 b 3 , y oc 0,3 + -\/3a4 + 63, with the usual normalization 

(78b 78+) = M 2 + M 2 + M 2 + H 2 + M 2 - 
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Tabic 1: SM singlet VEVs in our Higgs sector. 


label 

c PS 

C SU(5) 

C SO(10) 

c e 6 

label 

CPS C SU(5) C SO(10) 

c e 6 

Cl 

(1,1.1) 

1 

1 

27 

rfi 

(1,1,1) 

1 

1 

27 

c 2 

(1,2,4) 

1 

16 

27 

d 2 

(1,2,4) 

1 

16 

27 

ei 

(1,3,10) 

1 

126 

35T 

fi 

(1,3, TO) 

1 

126 

351' 

e 2 

(1,2,4) 

1 

16 

351' 

h 

(1,2,4) 

1 

16 

351' 

C 3 

(1,1,1) 

1 

1 

35T 

/3 

(1,1,1) 

1 

1 

351' 

e 4 

(1,1,1) 

24 

54 

35T 

h 

(1,1,1) 

24 

54 

351' 

es 

(1,2,4) 

24 

T44 

35T 

h 

(1,2,4) 

24 

144 

351' 

U\ 

(1,2,4) 

1 

16 

78 






u 2 

(1,2,4) 

1 

16 

78 






w 

(1,1,1) 

1 

1 

78 






V 

/ 

1 

45 

78 






y 

/ 

24 

45 

78 







(VEV terms only): 

J 35 i'x 35 V = 351'^a, 351T" = ei/i + e 2 / 2 + e 3 / 3 + e 4 / 4 + e 5 / 5 , (6) 

-^27x27 = 27^ 27 M = Cidi + C 2 d 2 , (7) 

^ 78 x 78 = 78^ 78"^ = 2mim 2 + w 2 + v 2 + y 2 , (8) 

7 3 5 i ' 3 = 351 ,/w * 351'^ 351 /A7 = 3 (e 3 e 4 2 + eie 5 2 - V2e 2 e 4 e 5 ) , (9) 

7^ = 35P Ma 35H^ 3M 7 A7 d a ^ #" A = 3 (/ 3 / 4 2 + fj 5 2 - V^/ 2 / 4 / 5 ) , (10) 

7 27 2 x 35 T 7 = 27 M 27" = c 2 2 / 4 + ^ 02/2 + c 4 2 / 3 , ( 11 ) 

J 27 2 X351' = 351'^" 27, 27, = d 2 2 e! + V^dirf 2 e 2 + d 4 2 e 3 , (12) 

hr 3 = 27^ 27" 27 a d^ = 0, (13) 

7 2?3 = 27^ 27j, 27 A cP" A = 0, (14) 

^27x78x27 = 27 M 78^ 27 

= Hg'UiCida + ^^ 202 ^! - (4cidi + c 2 d 2 ) - \ sj\vc 2 d 2) (15) 


7351'x78x 351' — 351' Mi , 78 K 351 7 — 

= ( V / 2 e 2 /i + V2e 3 f 2 + e 5 / 4 j + ^yg '«2 (v^ei f 2 + V2e 2 f 3 + e 4 /sj 

+ x^/ 2 u ' ( — 2e 4 /i — 5e 2 / 2 — 8e 3 / 3 + 4e 4 / 4 + e 5 / 5 ) — \\J\v (2e 4 / 4 + e 2 / 2 — e 5 / 5 ). 

(16) 

The general formula for ZT-terms in our case yields 

L> A = (27% (t A 27f + (27% (t A 27) M + (78 f )% (t A 78% 

+ (351%, (; t A 35T% + (35P%" (t A 351%,. (17) 
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Of the 78 D-terms, 5 vanish non-trivially, corresponding to the following generators of 
the SU(3)c x SU(3 )l x SU(3)# subgroup of E 6 : , t R , t R , t R , t R , given explicitly by 

D l = 73 0 Cl l J + l C2 l 2 + 2 l ei l 2 + 2 l e2 l 2 + 2 l e s! 2 _ I ^4 1 2 - |e 5 | 2 

_| dl | 2 _ |rf 2 | 2 _ 2\fi\ 2 - 2|/ 2 | 2 - 2|/ 3 | 2 + |/ 4 | 2 + |/ 5 | 2 ) , (18) 

-Dfj — | (—31 c 2 1 2 — 61 e 1 1“ — 3|e 2 |" + 3|es| 2 — |tti| 2 

+3|d 2 |“ + 6|/i| 2 + 3|/ 2 | 2 — 3|/ 5 | 2 + |w 2 | 2 ) , (19) 

D R — J 2 ^6c 2 ci* + 6cic 2 * — VSu\W* — VSwu7 + VSuiV* + VSvu\* 

—6d 2 d\* — Qd\d 2 * + v/3 u 2 w* + a/3 wu 2 * — V5u 2 v* — VSvu 2 * 

+6V2eie 2 * + 6^/2e 2 el* + Q\^2e 2 e7 + 6\^2e^e 2 * + 66465 * + 6 e 5 e 4 * 

- 6 ^/i/a* - 6 ^/ 2 // - 6 v^/ 2 / 3 * - 6 V 2 f 3 f 2 * - 6 / 4 / 5 * - 6 / 5 / 4 *) , ( 20 ) 

.Djj = ^ ^6 c 2 c 4 * — 6c 4 c 2 * — v/hiiin* + VSwui + VSuiV* — v^nni* 

+6d 2 ^i* — Qd\d 2 * — Vs u 2 w* + VSwu 2 * + V5u 2 v* — VSvu 2 * 

+6v / 2e 4 e 2 * — 6v / 2e 2 e 4 * + 6-\/2e 2 e 3 * — 6-\/2e 3 e 2 * + 6e 4 e 5 * — 6e 5 e 4 * 

+ 6 \/ 2 /i/ 2 * - 6 x/ 2 / 2 /i* + 6 \/ 2 f 2 f 3 * - 6 ^ 3 / 2 * + 6 / 4 / 5 * - 6/ 5 / 4 *) , (21) 

D R = 273 ( _2 l Cl l 2 + N 2 + 2 l ei l 2 _ I 62 ! 2 - 41 6 3 1 2 + 2 |e 4 1 2 - 1 65 1 2 + |'Ui| 2 

+2|di| 2 - |d 2 | 2 - 2|A| 2 + I/ 2 I 2 + 4|/ 3 | 2 - 2|/ 4 | 2 + |/ 5 | 2 - |u 2 | 2 ) • (22) 

They can be rewritten into 3 independent D-terms: 

D 1 = | Cl | 2 - |d 4 | 2 + | e 2 1 2 - I/ 2 I 2 + 21 e 3 1 2 - 2|/ 3 | 2 - |e 4 | 2 + |/ 4 | 2 - §M 2 + ||n 2 | 2 , (23) 
D n = |c 2 | 2 - |d 2 | 2 + | e 2 1 2 - I/ 2 I 2 + 21 e 4 1 2 - 2|/ 4 | 2 - |e 5 | 2 + |/ 5 | 2 - \\u 2 \ 2 + §|«i| 2 , (24) 

D 111 = +C 1 C 2 * — ^wui* + ^-vu\* + V2e 2 e \ 1 + V / 2e 3 e 2 * + e 5 e 4 * 

- d 2 d,* + f u 2 w * - f u 2 v* - V 2 / 1 / 2 * - V2/ 2 / 3 * - / 4 / 5 *, (25) 

via the definitions D 1 := VSD\ + 2Z)^, D n := —2D R and D 111 := D R + iD R , where D UI 
now forms a complex equation. The other independent combination D y = Z4| + a/3Z1|,+Z1jj 
is trivially zero, as it should be, since this O-terrn corresponds to the unbroken generator 
Of U(l)y. 

3.2 A specific vacuum solution 

In this section, we obtain a vacuum solution which breaks the gauge group to the Standard 
Model group SU(3)c x SU(2)^ x U(l)y. Due to the complexity of the equations we are 
unable to provide a full classification of vacua; a short discussion on alternative vacua can 
be found in Appendix |Dj 

The equations of motion in SUSY models are 

D a = 0, 

F<t> = 0, 
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(26) 

(27) 


with the usual definition of the F-term: 


FA := 


dW 


(28) 


In our case, there are 19 Standard Model singlets in the Higgs sector, giving 19 non¬ 
trivial F-terms, which can be easily reconstructed from the superpotential in equation (|5]) 
and the all-singlet terms of the invariants given by equations ([ 6 ]) (16). The non-trivial 
F-terms are given by equations (23)—(25). 


To obtain a vacuum solution, we perform the following steps: 

• First, we notice that F y leads directly to y = 0 (y is present only in the mass term 
78 2 ). By taking the self-consistent ansat^j] 


ci = di = f 5 = e 5 = 0 , 
U\ — U 2 — e 2 — f 2 ~ 0, 


(29) 

(30) 


the system of equations is greatly simplified: F C1 j Fd i) F e2 , Ff 2 , F e5 , Ff 5 , F Ul , F U2 , F y 
and D ni are solved automatically. 


Solve F ea and Ff 3 for f 3 and e 3 , respectively to get 

9A 2 / 4 2 


e 3 = 


3m 35 r - y/2\%w' 


fz — 


9Aie 4 2 


3m 35 i' - V2\$w 


(31) 


Solve F C2 and Fd 2 for fi and e 4 , respectively to get 

C 2 ( fK\ ( fTF.. , . 


ei = 

fi = 


24A 4 d 2 

d 2 


- 12777.27J , 

(32) 

- 12 m 27 ^ . 

(33) 


24A 3 c 2 

Simultaneously solve F ei and F f 4 for U 

_ (3\/2m 35 i/ - 2A 8 w)(3v / 2m 35 i/ + A 8 w) 

324e 4 A 4 A 2 

Simultaneously solve F ei and Ff t for d 2 \ 

(6y / 2/n 2 7 — A 7 (V / l5^ + uj)) ( 6 ^^ 351 / — As(\/l5u + ta)) 


(34) 


do — 


144A 3 A 4 c 2 


(35) 


It is now convenient to define a new quantity A := \/l5u + w. We can now solve F v 
for v as a linear equation: 


v = 


,MA 2 \$\s — 8v^5A7j 4(??7. 3 51' A 7 + 2m 2 7As) + 24\/l0m 2 7(2m 3 5i/A7 + m 2 7 As) 


768-\Z3ru7 8 A 3 A 4 


( 36 ) 


^^This ansatz is motivated by the ansatz used in the model with the 78 omitted m , where a classification 
of vacua is known. See Appendix |D| for further discussion. 














• Three variables remain to be determined: A, c 2 and e 4 . We are left with only one 
unsolved F-terrn F w , which is a polynomial in A: 

0 = P 0 + Pi A + P 2 A 2 + P 3 A 3 + P a A\ (37) 

Pq — — 576 m27 (2777357/A 7 + 777.27 Ag) ^ 25 m 2 Y (2777357/A7 + 77727 Ag) Ag 

— 480?77.35^/m7gA3A4Ag T 110592??r^gAiA2A3A4^, (38) 

P 1 = 192V2(995328X 1 X 2 XlXlm 3 78 

+ 46O8A3A4 ( 24 m 27 AiA 2 A 7 A 8 + m354/ (12A1A2A 2 — AgAiAg)) r?7 2 g 


+ 240 A 3 A 4 Ag (—m3 51 /Ay + 2777357 / 7772 7 AgA 7 + 2 ?? 7 2 7 Ag) 7 ? 7 7 g 
+ 2577 l 2 7 A 7 Ag (2777.354/ A7 + 5777357/77727 AgA 7 + 2r?72 7 Ag) ^, ( 39 ) 

P2 — —16 ^Ag (I8432A3A4 (9A1A2A7 + A 3 A 4 Ag) 777^3 
+ 240 A 3 A 4 A 7 Ag (5777351/A 7 + 16 m 27 A 8 ) 777 78 

+ 25 A 7 Ag (2777354/A7 + 14777354/77727A 8 A7 + llm^Ag)) ^, ( 40 ) 

P3 — 40 v / 2 A 7 Ag ( 96 ?t 7 78 A 3 A 4 + 5A7 (m35i/A7 + 2m 2 7X %)), ( 41 ) 

P 4 = — 25 A 7 Ag. ( 42 ) 

Note that the coefficients Pj depend only on the Lagrangian parameters; choosing 
those, we can determine A numerically. 

Finally, we solve the remaining D-terms and determine C2 and e4 from 

0 = D 1 = 2 |e 3 | 2 - |e 4 | 2 - 2|/ 3 | 2 + |/ 4 | 2 , ( 43 ) 

0 = D n = |c 2 | 2 - |d 2 | 2 + 2 |ei | 2 - 2 |/i| 2 . ( 44 ) 

Note that v and w are determined, once A is fixed. We can therefore see that 


f A oc l/e 4 , e 3 oc /| oc l/e| and f 3 oc e 2 . Therefore D 1 can be written as a quartic 
polynomial in |e 41 2 ; the constant term has a positive coefficient (the 63 term), while 
the highest order term in |£ 41 2 has a negative coefficient, a solution for a real 64 > 0 
will always exist. Similarly, d 2 oc l/c 2 , e 4 oc c 2 /d 2 oc c| and f\ oc d 2 /c 2 oc l/c|, D n is 
a quartic polynomial in C 2 independent of e A \ the constant coefficient will be negative 
(the f\ term), while the highest order coefficient in |c 21 2 is positive (the e\ term), 
which again guarantees a real solution C 2 > 0. 

The initial Lagrangian parameters constitute of the 3 masses 77727 , 777357 /, m 78 and 8 
massless parameters A,, i = 1,... , 8 . The simplest order in which to compute the given 
vacuum solution with these parameters is given below: 

1 . Take the ansatz for some of the VEVs (symmetric under conjugation symmetry [IT] ) 

Ci = 0, d\ — 0, 

e 2 = 0, f 2 — 0, 

65 = 0 , /s — 0 , 

771 = 0 , u 2 = 0 , 

y = 0 . ( 45 ) 
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2. A is determined through the polynomial, and then the VEVs v and w are determined 
by 


v = 


x/X0H 2 A 2 A 8 — 8\/5A7^4 (77T351'A7 + 2 rri 2 T \%) + 24\/l077727 (2777351/A7 + ?7727A 8 ) 


w = A + 


768^^78^3-^4 

— 24*^/277727 (2/77351' A7 + 77127X3) 


76877178X3X4 

■h A7-A (8777351'A7 + 1677727-^8 — V / 2AA 7 A 8 ) 


(46) 


(47) 


3. 

4. 


e 4 and C 2 are determined through D 1 and D n , respectively, using equations (48)—(53). 

The remaining nonvanishing VEVs are o^, e 4 , fa,e 3 , fa, fa, and they can be computed 
in terms of A, w, C 2 , e 4 : 


(48)—(53 


C?2 


ei 


fa 


e 3 


fa 

fa 


144A 1 3 A 4 c 2 ( 6v ^ m27 “ A? ' 4 ) ( 6v ^ m351 ' “ Asj4 ) ’ 
I 2 A 3 C 2 2 

a/2AA 8 — 12/77351' 

^AgA 1 !” 27 ( 6 ^ m27 " MA ) ( 6 ^ m351 '" As / ■ 
(3^777351' - 2A 8 U/) (3V^r77 35 i' + A 8 w) 


5832 v / 2e 4 2 AfA 2 


9e 4 2 Ai 


3777351' — 

1 

324e 4 AiA 2 


\J~2w A 8 ’ 

( 3^777351' - 2A 8 7 V ) ( 3-^/2777351' + A 8 W7 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 


The solution above does indeed break E 6 into the Standard Model group. This can be 
checked by explicitly computing the gauge boson masses, found in Table [5] of Appendix |Aj 
It is possible to further illuminate this result by considering that under the standard em¬ 
beddings of the E 6 subgroups from Table [lj the SU(5) breaking in our solution is solely 
due to non-vanishing VEVs e 4 and fa. 


4 Doublet-triplet splitting 

We tackle now the issue of doublet-triplet splitting. We denote the doublets and anti¬ 
doublets by D ~ (1,2,+1/2) and D ~ (1,2,—1/2), while the triplets and antitriplets are 
denoted by T ~ (3,1, —1/3) and T ~ (3,1, +1/3). The detailed labels of these states are 
given in Table [6] of Appendix [B] 

I 11 the Higgs sector of our model, there are 12 doublet-antidoublet pairs and 13 triplet- 
antitriplet pairs. This is one extra pair of each compared to the renormalizable model with 
the Higgs sector 351' + 351' + 27 + 27, in which doublet-triplet splitting surprisingly fails 
(see m for details). The extra states come from the added 78, and are labeled by the 
index 0, i.e. Dq,Dq,Tq,Tq. We shall see that this extra row and column, together with a 
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new vacuum compared to the model without 78 now enable doublet-triplet splitting in the 
usual way (by fine-tuning). Note that all the doublets and triplets are located in a 5 or 45 
(or their conjugates) of SU(5), except for one extra triplet in the 50 of SU(5). 

The mass terms for the doublets and triplets are written as 


Do ■ ■ ■ Dii ) M 


11 I •'' /t doublets 


/DA 


\DnJ 


+ Tn ■■■ TAM 


12 I triplets 


/TA 


\tJ 


(54) 


The matrices A4d ou biets and A4 triplets are similar; we can compactly write a 13 x 13 
matrix A4 with block form 


M = 


Mu Tf 1 2 

M‘2l M‘22 


(55) 


with the diagonal blocks defined by 


Mu = 


M‘22 — 


/ 

77178 

0 

2%/3 

0 

0 


0 ^ 




0 

m27+M ^V5, 

3/4 

vT5 

6 A 5 C 2 

0 


0 




2v/3 

aA 4 e 4 

s/15 

m27+M ^M 

0 


0 




0 

6Agd2 

0 

. A 7 (3v'3v — 
m27 + 1 6/15 

w) 0 


_ vTo 

5 



0 

0 

0 

0 

, A s n A S UI 
m 351 ' + 51/M S/I 


-y/|a!Aie 4 



V 

0 

0 

0 

_ a 3 c 2 

\/Io 

^|aA 2 /4 

rr, A 8 t’ >~8 w 

m 35i' S/W S7I / 



/ 

m 351 

-It/IO 8 12V2 

0 

0 

|a/3Q!A2/4 

0 


0 \ 



0 vspMfi 

0 

0 


0 

0 



0 

\/5/3Aie 4 

m 35T + i^ _ ^i 

0 

0 


0 

2\/T0Aie 4 



3o:Aie 4 

0 

0 m 351 , 


0 


0 

0 



0 

0 

0 

0 m 351 / 

4v / 30 + 12^2 

0 

0 

V 


|\/l5/3Aie 4 

0 

0 

0 

0 

’"351'-4^B A 3 , ' + 1S/3 

0 


0 

0 

SV55A 2 /4 

0 

0 


0 

m 35T - TAh - 677f / 


(56) 


and the off-diagonal blocks defined by 

/ * 8 /3 

2 a/6 

M\ 2 = 


, (57) 


M 2 i — 


2v/6 

0 

0 

0 

0 

0 

A 8 e 3 

2\/6 

0 

0 


0 


A S e 4 _bsJ± VSp\)je 4 q \ 

!4a/2 2a/6 24a/2 


24a/2 

0 2 V ^A 3 C2 
0 0 

0 0 


A S«1 
2\/6 

24a/2 

0 0 


2\/6 

0 

-\/2A 4 d2 

0 

0 

0 


24 a/2 

0 

0 

0 

0 

0 


0 

0 

0 

0 / 

o/ 


(58) 


0 0 0 \ 
— \f\^302 0 0 


0 


0 


-\/2A 3 C2 


0 0 
0 0 

0 0 

0 0 
0 0 . 


(59) 


The separation of the matrix A4 into 4 blocks is arbitrary and is used above merely 
as a simple way of presenting a large matrix. The matrix A1 doublets is obtained out of A4 
by removing the last row and column, and taking a = —3, (3 = —a/3, while the matrix 
A4triplets is obtained by taking a = (3 = 2. Note that the matrices were already simplified 
by taking the vacuum ansatz of vanishing VEVs 


Ci — di — e 2 — f 2 — 65 — /s — Ui — U 2 — y — 0, 


(60) 
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while the rest of the vacuum solution was not plugged-in. Notice also that the coefficients 
a and (3 are located as factors in front of e 4 or / 4 ; this is expected, since e 4 and / 4 are 
the SU (5)-breaking VEVs, so they control the difference between the doublets and triplets. 
Also, the coefficients a are —3 and 2 for the doublets, respectively, which come from the 
VEV (24) in the terms 5 • (24) • 5 of SU(5). The coefficients /3 are — \/3 and 2 for the 
doublets and triplets respectively, which are the Clebsch-Gordan coefficients coming from 
the terms 5 • (24) • 45 or 5 • (24) • 45 of SU(5). One can check that this is indeed the case by 
locating the position of a and (3 coefficients in M. and checking, to which states in Table [6] 
these matrix entries correspond to. 

The procedure for fine-tuning is now similar to the one attempted in the 351/ + 351' + 
27 + 27 case in |jTT] . Once the vacuum solution is plugged-in, the determinants of the two 
matrices become zero: 


det ( Af doublets ) 


vacuum 


0 . 


(61) 


det (Aftriplets) 


vacuum 


0 . 


(62) 


The massless doublct-antidoublet and triplet-antitriplet pairs are simply the unphysical 
would-be Goldstone bosons from the breaking of Eg. They correspond to the broken gen¬ 
erators with the same quantum numbers, which can be found in SO(IO) language in the 
16 and 16 parts of the adjoint 78. The doublet-triplet splitting condition then requires 
another doublet-antidoublet pair to be massless, which imposes the following condition: 


Cond(A4) := pmin, M ., 


(63) 


where pmin,- denotes the i- th principal minor of rank 1 (the sub determinant when the i-th 
row and column are removed). In our specific case, we chose to remove the first row and 
column, which correspond to the fields D 0 and D 0 . We chose these due to simplicity, but 
note that this choice is valid only if the Goldstone modes of the given vacuum have nonzero 
Do and D o components. This can be checked later on via equations (71) and (72). The 
logic behind the whole method is explained in Appendix [Cj 

Given the notation above, DT-splitting can be performed by a fine-tuning, so that 


Coild (-Ad doublets) 0, 
Cond(A4 triplets) ~f~ 0. 


(64) 

(65) 


The above conditions can in principle be computed analytically, but they are too compli¬ 
cated to be of any practical use. The viability of the splitting conditions can be shown 
numerically though, as well as with some careful considerations. The most convenient way 
to perform the DT splitting is to make use of the dimensionless parameters A 5 and Ag from 
the superpotential (see ([ 5 ]) ). Since the invariants 27 3 and 27 have no all-singlet terms, 
A 5 and Ag are not involved in the computation of the vacuum. Moreover, the matrix A4 
contains only a single entry with the parameter A 5 and a single entry with the parameter 
Ag. The condition in equations (64) and (65) can be written as 


K\ — K -2 A5 Ag — 0 , 
K[-K ' 2 AsAg^O, 


( 66 ) 

(67) 
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where K±, K 2 , K [, AT, depend only on the other parameters in the Lagrangian (771351/, 77773, 
77727, Ai, A2, A3, A4, A7, As) and the vacuum itself, all of which can be computed inde¬ 
pendently from the parameters A 5 and \q. The analytic forms of K s are not very illu¬ 
minating, but they can easily be evaluated numerically for any values of the parameters 
777351/, 777 78 , m 27 , Ai, A 2 , A 3 , A 4 , A7, Ag. One can then take 


A5 — 


K\ 


K 2 \ & 


( 68 ) 


with Ag arbitrary, and get an extra massless doublet mode. Substituting into condition (67), 
we can numerically check that indeed 


K[K 2 ± K’ 2 K x . 


(69) 


We therefore conclude that the addition of the 78 in the Higgs sector, both with the 
new vacuum and enlarged matrices Afdoubiets an d AT triplets, now allows for a DT splitting 
in the model. Without the 78, a similar attempt at fine-tuning is not possible, since in 
that case we get K 2 = K' 2 = 0 after inserting the vacuum and fine-tuning in A 5 or Ag is not 
possible. Note that this tree-level fine-tuning is stable under quantum corrections due to 
the non-renormalization theorem for the superpotential. 

The linear combinations of Ds and Ds , which correspond to the MSSM Higgses H u and 
Hd, can be found by computing the new left- and right- null-eigenvectors of the fine-tuned 
matrix ATdoublets- In practice, the physical MSSM Higgses can most easily be extracted by 
computing the (left and right) null-eigenspace of the fine-tuned matrix AT doublets, which is 
2-dimensional. Given any basis of the null-eigenspace, the Higgs will always be orthogonal 
to the would-be Goldstone boson. The would-be Goldstone itself can be easily identified 
by the fact that it has no component in the directions Di (or Dj) for i = 1,3,4, 5, 7, 8 . 
This absence of some doublets in the Goldstone can be deduced from the mass matrix 


in equation (55), but we also confirmed this by noting that the Goldstone components (f>i 


are the ones which have couplings of their derivatives to the the gauge field through the 
following type of expression (originating from the kinetic terms of scalar fields): 


-igK W ) (*“W)‘. 


(70) 


Choosing a to be the doublet/antidoublet broken generators and using our vacuum, we 
identify the doublet components in (j)\ to which there is no coupling as indeed those listed 
above. Explicit computation identifies that the prevailing cause of some components not 
being present in the would-be Goldstone mode is our ansatz of vanishing VEVs. More 


precisely, without the vacuum ansatz, the terms in expression (70) can be schematically 
written as 
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dDh k + 3V5v + 5V3w-Vm, gD , _d laD ,^d 1 gD , + e, /3 ^ 


<+D G * oc 


10 


_ 7Y c®; _ £l an; _ 7+ d o' 7 _ ll.fi go; 

4 5 6 4 7 4 V 2 8 


10 a/2 

i V^5/2 an* | /3 qn* i V^8,/2 o n* . ^5 /h an* 

4 a/2 4 7 4 V 2 8 

- e i^ dD ° + T2 dDl °- e i^ dD "' 

3y/5v + 5a/3w — a/30 y —* c 2 — * C\ n -=r* /s /3 , — 


<+D n + -f dD, + -+<TD, ~ H-\ -dD 


( 71 ) 


(72) 


We see that the absence of components Di and D, t for i = 3,4, 5, 7, 8 are directly related to 
the vacuum ansatz with vanishing VEVs. D\ and D\, however, are always absent])] 

The ratios of the various components of the Higgs, and hence the ratios of the EW 
VEVs of these component, are computed from the null-eigenspaces, while the magnitudes 
of the VEVs are given by the VEVs v u = (H u ) and Vd = (H d ). We have 


v l = 

(73) 

7=0 


11 


v d = X]l+I 2 > 

7 = 0 

(74) 


where ty = (Di) and Vi = (Di), and additionally the following MSSM relations hold: 

(75) 

(76) 


v u + v d ~ (246 GeV) 2 , 
Vu/vi = tan 0 . 


5 Yukawa sector 

The Yukawa sector comes from the Yukawa part of the superpotential after inserting the 
vacuum solution: 


IVYukawa = 27/ 27/ (YU (27) + Yjp (351')). (77) 

In addition to the GUT scale VEVs in the representations 351' and 27, the EW Higgses 
H u and Hd also need to be present in both of these representations. The low-energy MSSM 
Higgses H u and II,/ come from a linear combination of the doublets of type D r\_/ (1,2,+1/2) 

1 The 27 contains 2 singlets, but 3 doublets, so the singlets cannot be pushed to all the doublets with 
a t° generator, and thus one doublet component is missing. The representation 351' has 5 singlets and 8 
doublets, but the projection relation dijk 351° fc = 0 reshuffles the definitions of singlets and doublets, so 
there is no loss of doublet components. 
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and antidoublets of the type D ~ (1,2, —1/2), respectively. These states are identified in 
Table[6j and their EW VEVs are labeled by ty := (Dj) and bj := ( D t ), where i — 0 , 1,..., 11. 
As discussed in section |4j we generically expect ty ^ 0 and v t ^ 0. 

Each generation of fermions is found in a fundamental representation 27 of Eg. The 
matter content of each generation, subdivided into SO(IO) representations, is the following: 

• The 16 of SO(IO) contains the SM particle^/] and the right-handed neutrino u c . 

• The 10 of SO(IO) contains a vector-like pair of down-type quarks (3,1,—1/3) + 
(3,1, +1/3), as well as a vector-like pair of lepton doublets (1, 2, —1/2) + (1, 2,1/2). 
These exotic are labeled d', d ,r \ IJ and L" , respectively. 


• The SO(IO) singlet 1 is denoted by s is sterile, and has a role analogous to the 
right-handed neutrino. 

Suppressing generation indices, the mass terms are computed to be 

U T (-Ul+27 + 

+ (d cT d ,cT ) ( l * * ’ 2Y27 + +7To V4 ^51' + 276^35+ C2 * 27 + 775^35F 

V ) \-V 3 Y 27 - ^VgY^r ~ ^VU -Cll'27 + 

-V 2 Y 27 - ^—74^3517 + ^/|ws%F C 2^27 ~ |^h/5%F 

27 + ^TTo^^lsl 7 — ~ Cl ^27 — 2 7TH^ 4 ^l5T 

t ’ 1 ^ 27 ~~ 27To W 5 ^35l 7 — 

_ 72 t ’ 1 °-^35T 7 -V1Y27 

1 /1%F ^hY. 

,cT ' y / 2 % f / 3 %f S2V27 - 


+ (e T e /T 


+ ( Z/ T l/ T 



8 V7 ^351' 


'v^ 6 ^ 351 ' 


C2Y 2 7 - 


2 +I5 ji> - 1 351' 

3 1 


5 V 5^35l' C1I27 2 4 ^351 


- I 7/ cT S 7 I/" 


—V3Y27 + \/lv 9 yw?\ 


/« 


77 


s 

/C 


351 

-V 3 Y 27 + x/lbgy+iT U2+27 


V I Z/ T i7 ,T 


V 

^1 Y 35V 

,72^2^35+ A 3^%F 


351' W£l \l 5 U 4 E 35 I 7 


351' 

2 - 


0 


7 


/ Z7 C \ 

s 

w 

(78) 


Notice the different Clebsch-Gordan coefficients in front of iy (i = 2,3,4, 9) that come 
from the couplings with the 5 of SU(5) and in front of bg^i that come from the 45, while 
C 12 originate from a SU(5) singlet and from a SU(5) adjoint 24. 

The A i correspond to VEVs induced in the electrically neutral component of weak 
triplets of type (1, 3, +1), while Aj correspond to VEVs induced in weak triplets (1, 3, — 1). 
The Aj and Aj are found only in the representations 351' and 351', as shown in the defi¬ 
nitions of Table [7j The mass matrix has contributions from the terms 351' • 351' and 

1 We use the standard notation Q , L, u c , d c and e c for the representations (3, 2, +1/6), (1,2,—1/2), 

(3,1,—2/3), (3,1,1/3) and (1,1,1), respectively. The lepton doublet contains the electron e and the 

neutrino ta 
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351' • (78) • 351'. Its explicit form is computed to be 


M a = 


^m 3s v 

(A 2\/A + 5 vl) 

ujAg 

2\/3 

0 

6Aie! 

\ 


2^“ 2 

m 351 ' + A 8 ( 12v /2 4 a /30 _ 2 vl) 


—6Aie2 



0 


m 3 si' + A s (545 - ^ - |fj) 

6Aie 3 

17 s) ) 

V 

6A2/1 

~ 6A2/2 

6A2/3 

TO351' + As (573 + tIj + 


( 79 ) 


Integrating out the heavy weak triplets and inserting the F-term ansatz of vanishing VEVs, 
we get 


/AA 

A 2 

A 3 

w 


As (e !/2 2 

0 
0 

6 A 2/1 


0 0 

d +A8 (l^f + i^s) 0 

0 m 35i'+ A s(^-^) 

0 6 A 2/3 


6Aiei 

0 

6 Aie3 

' ,+As (371 


/ A 4 u 3 2 \ 
X 4 \Z 2 v 2 v 3 
A 4 u 2 2 
A 3 V 


(80) 


After integrating out the heavy vector-like states from equation (78), and using the 
ansatz c 4 = f 5 = 0, we get the matrices for the low energy states: 


M v = - Vl Y 27 + (^u 5 - ^u 7 ) Aoio (81) 

(1 + (9/4)XXt)" 1/2 ((a 2 - lv 3 X) Y 27 + (^4 - \v*X) + ^(t/ 8 - ftJnX)) Y^r), 


M t d = 


m e = { 1 + xxt) 1/2 ( ( —,- 2 - v 3 x) y 27 + + v 9 x) + ^(vs + €nX)) i^), 

m n = — (i + xx*y 1/2 


1 Vltlfi / 3 VIV7 I 1 UBWIO C2 1 . /c 4>7t/in C2 1 4 V 5 V 6 C2 O / 10 A C2 \ V 

_ VnTT ~ V 2T +75TTJI + vs TrE + TOTS ~ WT A2 /7 I 127 


(82) 

(83) 


2 /i 73/1/4 


h h Ys h h 


h 


3 U 7 U 5 3 n 7 2 1 n 6 2 

+ + 7T—- A 1 I 135U 


1 V5 2 

40 /, V 80 /, 8 /, 2 ,/ 3 


I I vY_ _ 9 IOV 1 V 1092 10 T'lQ 2 C 2 2 I 9 /10 1>1 U 6 C 2 I 8 /F £ 2 l _ 20 A \ y V ' -1 V 

+ ^ fi 2 V 3 h U + 3 fi fY + Z \f 3 f 3 h + 3 /a / 4 2 3 As / 4 2 J ^Asl^ 27 


/8\/l0 


Ult> 5 C 2 


V 3 h fY 


Y 27 Y—Y 27 Y—Y 27 

27 351' 27 351' 27 


-- A ) Y 2T Y^,Y„Y^,Y„Y^,Y„ 


3 / 3 A 2 
(i + x*x T )~ 1/2 


351' 


where 


x = -^IItY 27 y^± 


h 


(84) 


(85) 


Notice that the main factor in the expressions is a linear combination of the following 
matrices: 


Ur, %17, Y 27 Y^LY 27 , Y 27 YY±Y 27 YY±Y 27 YY±Y 27 , (86) 


' 351' 


351 
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which are all symmetric, since Y 2 7 and are symmetric, as are Mjj and Mjy. The 
matrices Me and Me are not symmetric, though, due to the projection factor onto the 
light families. 

As usual with vector-like states, the expressions for the low energy masses in equa¬ 


tions (81)—(84) are nonlinear, which complicates the analysis of the masses. We comment 
more on the low-energy part of the Yukawa sector in section [ 6 j in which we also do a 
numeric fit in the 2 -generation case. 


6 Numeric fit of the Yukawa sector 

The presented model has 3 masses and 8 dimensionless parameters A* in the breaking 
sector, as well as two symmetric 3x3 Yukawa matrices in the Yukawa sector, which is 
easily seen from the superpotential in equation ([!]). Note also that the product A 5 A 6 is 
fixed by the fine-tuning of the EW Higgs mass. Also, a rotation in family-space can bring 
one of the Yukawa matrices to be diagonal. We shall limit ourselves to the case where all 
parameters of the Lagrangian are real. The independent number of real parameters in the 
breaking sector is thus 3 + 8 — 1 = 10, while the 3-family Yukawa sector has 6 + 3 = 9 real 
parameters. Since the number of independent real parameters is 19, while there are only 
17 real numbers to be measured (3 masses in the up, down, charged lepton sector each, the 
two differences of masses-squared in the neutrino sector, as well as 3 angles in each of the 
CKM and PMNS matrices; we neglect the CP-phases), the general expectation is that a 
fit is possible to perform. There may exist, however, non-obvious mass relations concealed 


due to the complexity of the low energy expressions in equations (81)—(84), which are not 


respected by the experimental values; these fears can be alleviated by finding points in 
parameter space, which give a good fit for to the masses and mixing angles. 

In this section, we perform a fit in the simplified case of 2 families. Here, there are 
again 10 real parameters in the breaking sector, while the Yukawa sector has 3 + 2 = 5 real 
parameters. The fit is performed for the results of the quark masses m t , m c , mb, m s , the 
charged lepton masses m T and m M , the difference of the squared neutrino masses m^ 3 — r ml 2 , 
and the mixing angles 9 c b and 6*23 in the CKM and PMNS matrices, respectively. 

The most convenient way to perform the fit is to take some of the GUT scale VEVs 
as the parameters in the fit, instead of the initial parameters in the Lagrangian. This 
is advantageous since the equations of motion are linear in the Lagrangian parameters. 
Taking the ansatz <+ = cl\ = e 2 = fi = e 5 = / 5 = u\ = u 2 = y = 0, we need 12 parameters; 
only 10 are independent, while 2 are determined through the Z+terms. It is convenient, for 
example, to take the following quantities as independent parameters: 


^351', C 2 , d 2 , /l, / 3 , fi, e 4 , V, w, A 


5- 


(87) 


We can then use the equations of motion to determine the remaining initial parameters 
and VEVs. The ZA-terms are solved by taking 


ei = ± v + 2 | 2 - c 2 | 2 + 2|/i| 2 /V2, 

<•3 = ±VM 2 - l/il 2 + 2IAI7V2. 


( 88 ) 

(89) 


Since we want ei and e 3 to be real numbers, the arguments in the square root need to be 


positive, which limits the space of parameters in equation (87). Alternatively, we could also 
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take the 6 independent parameters, for example^ to be the VEVs d 2 , fi, e 3 , / 4 (now with 
no restriction) and two angles 61,62 G (0,27r), with the remaining quantities computed as 


fi = sin &i a/ c 2 

2 + 2 e 4 

vA 

d 2 = cos 61 ■>/ c 2 

2 + 2 e 4 

12 

1 ? 

/ 3 = sin do a/ / 4 

2 + 2 e 3 

7 A, 

e 4 = cos 6 2 Vl /4 

2 + 2 e 3 

2 


The F-terms then yield 


(90) 

(91) 

(92) 

(93) 


m 27 = 


^351' 


2 v(w + VTEv) (e 4 / 4 - 2 e 3 / 3 ) : 


m 78 

Ai 

A 2 

^3 

a 4 

A7 

As 


2 w (a/15w - n) c 2 d 2 (4 e 3 / 3 + e 4 / 4 ) 

+ e i/i(' y — T/lSta) ^2e 3 / 3 (3A/l5n — 5w) — e 4 / 4 ( 3 a/15u + 7 m;) 

\Zl5m 35 i/ (e 4 / 4 - 2e 3 / 3 ) 1 2 
2w (a/15w - w) (4e 3 / 3 + e 4 / 4 ) ’ 

/ 3 /4^ 3 51' 

e 4 (e 4 / 4 + 4e 3 / 3 ) ’ 
e 3 e 4 m 3 5i/ 


/ 4 (e 4 / 4 + 4e 3 / 3 ) ’ 

eim 35 i/ (a/X5u (e 4 / 4 - 2 e 3 / 3 ) + 3w ( 2 e 3 / 3 + e 4 / 4 )) 

2wc 2 2 (4e 3 / 3 + e 4 / 4 ) 

/i?n 35 i/ (\/l5u (e 4 / 4 - 2e 3 / 3 ) + 3w (2 e 3 / 3 + e 4 / 4 )) 

( 2wd 2 2 (4e 3 / 3 + e 4 / 4 ) 

3 ^/ 77351 / (2 e 3 / 3 - e 4 / 4 ) ((m - a/15w) e 4 A + 2v (2 e 3 / 3 - e 4 / 4 )) 
w (Vl5 w - u) c 2 d 2 (4e 3 / 3 + e 4 / 4 ) 

3 y / 2 ? 77 . 35 i' (2e 3 / 3 - e 4 / 4 ) 
w (e 4 /4 + 4e 3 / 3 ) 


(94) 

(95) 

(96) 

(97) 

(98) 

(99) 
( 100 ) 
( 101 ) 


Using these values, Ag can then be determined by hne-tuning in the doublet mass matrix 
(see section [4]), from which also the EW VEVs u,- and hj are computed. One can then use 
these to compute the mass matrices from equations (81)—(84), and ultimately the masses 
and the mixing angles. For our numeric fit, we assume: 


• the simplified case of second and third generation only; 

• all parameters real; 

• all errors in measured quantities at the 10% level; we believe that at the present stage 
such a choice is a good compromise between the realistic case and the simplicity of 
the analysis; 

• the values for the masses and mixing angles at the GUT scale ~ 10 16 GeV as shown 
in Table [d] [26] . valid for tan/3 = 10; 

1 In both parametrizations of encoding the F-terms, we wrote them so that we retained control to make 

fi and / 3 potentially small; this will be important for the neutrino sector, as described later in this section. 
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positive signs in eqs. (88)-(89). 


We performed the fit by minimizing the chi-squared function 



(fi(x) ~ Vi) 

4- of 


( 102 ) 


with respect to the initial parameters x. The experimental values are denoted by the 
values computed using our model are fi(x), and <7j are the 1-sigma deviations from the 
values y % . For our fit, we have 9 different measured quantities, thus i = 1,..., 9. The errors 
are taken to be cq = 0.1 y i: giving 

fj{x) -Vi V 

Vi J 

We define the convenient measure y 2 = y 2 /9, which tells us the average a 2 deviation per 
measured value. Also, we define the pulls Xi — (/*( x ) — 2 /i)/, which tell us how many 
sigma a certain quantity deviates from the measured one. 

Due to the large number of parameters, the y 2 function will have many local minima in 
the parameter space. We give below two such points in the parameter space, corresponding 
to the best fits that were found and which we deem sufficiently good (with y 2 < 1). The 
results are given in Tables [2| [3j [4j the points in the parameter space are given in terms of 
the independent parameters best suited to a numeric search (as discussed in this section) 
in Table [2] and in terms of the original Lagrangian parameters in Table [3j while Table [| 
shows the obtained results for the masses and mixing angles. 

Notice from Tables [2] and [3] that of the two Yukawa matrices Y 2 7 and Yyy, Yyy was 
chosen to be the diagonal one. Furthermore, the original Lagrangian parameters in Table [3] 
are given so that it can be checked they roughly fall into the perturbative regime. A possible 
exception could be the value As of the first solution. 

Note that the given two points are merely the best ones we found. Due to the high 
dimensionality of the parameter space, we suspect there are likely many more points which 
give a comparable or a better fit. There are a number of observations that can be made 
about these points in general, however, by deducing them from the formulae, comparing the 
two parameter points given in the tables and through experience obtained by performing 
the fit: 



(103) 


• The parameter points are not necessarily close to each other, which also holds true 
for any specific single parameter. We can see in Table [2] that it is not necessary for 
any parameter to be at a very specific value to obtain a good fit. The suitable areas 
of parameter space thus form many disconnected regions, and no specific value can 
be claimed for any parameter. In this sense, the mass formulae of the theory are not 
very predictive of the original parameters. 

• There is no specific mass or mixing angle, where one would consistently be getting 
tension. As observed from Table [4j while a specific solution might have most tension 
with observation coming from a single mass or mixing angle, there are other points, 
which also give a good fit and where this parameter is predicted better. In this sense 
we cannot claim any tendencies in tensions of the observables. 
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Imagine that in Table [2] or [3] we rescale all the mass parameters with a common factor. 
Such a rescaling would have no influence on the masses of the quarks and charged 
leptons, which are controlled only by the EW VEVs, as can be seen from equa¬ 
tions (81)—(83) and equation (85). Rescaling would influence the neutrino masses, 
however, due to the seesaw mechanism, confirmed by equation (84). In principle, the 
rescaling factor can always be adjusted, so that the fit of the mass-square difference of 
the neutrino masses ml — ml 2 is exact (provided this does not spoil the GUT scale, 
or the upper bounds on neutrino masses). Thus only neutrinos are actually sensitive 
to the GUT scale. 


• The mass parameters are chosen to be at approximately the GUT scale, say at the 
order of 10 15-16 GeV. If this is true for all the mass parameters, there might be a 
problem with the neutrino masses. It is a well known fact that this GUT scale is 
a few orders of magnitude too large compared to the seesaw scale for sufficiently 
large mass differences in the neutrino sector. One can cure this problem by having 
the spontaneous symmetry breaking occur in multiple stages, with a mass hierarchy 
between different VEVs and the seesaw scale corresponding to one of the intermediate 
stages. The gauge coupling unification in such a scenario could then be spoilt by the 
particles appearing at these intermediate mass scales. 


In our case, however, the number of different VEVs is large enough, so that having 
one or two of the VEVs at a smaller scale does not disturb the breaking pattern. 
To see this, note that to obtain sufficiently large neutrino masses, it is enough for 
one of the terms in equation (84) to be of the proper scale, which can be achieved 
by simply taking f\ or / 3 to be several orders of magnitude smaller than the GUT 
scale (the seesaw type I contributions). In the solutions given in Table [2j we achieved 
sufficiently high neutrino masses by taking the parameter f\ (and also v) to be a few 
orders of magnitude below the GUT scale. As seen from the gauge boson masses 
in Table ([5]), small f\ and v do not spoil the one-stage breaking scenario; since the 
intermediate scales are not associated to an intermediate symmetry breaking at that 
scale, but are instead purely accidental due to a carefully chosen parameter point, we 
do not expect too many relevant states (apart from the singlets u c with f\ Majorana 
mass) to be found at the scale /j. We noticed however the appearance of a lighter 
color triplet-antitriplet pair, possibly due to the similarities between the doublet and 
triplet mass matrices. Although this influences the running of the gauge couplings, 
we will neglect it in view of the (presumably) large threshold uncertainties present 
anyway. 


• A final comment on the neutrino masses: although only the difference of masses- 
squared needs to be fitted, one still needs to check that the neutrino masses themselves 
are < 1 eV m- As seen in Table [4j this condition holds true for both of our parameter 
points. 
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Table 2: Two example points in the parameter space, written in terms of the independent 
parameters suitable for a search. 


parameter point 1 point 2 


^ 35 i' [GeV] 
c 2 [GeV] 
d 2 [GeV] 
fi [GeV] 
h [GeV] 

U [GeV] 
e 4 [GeV] 
v [GeV] 
w [GeV] 

As 

(^27)11 
(^27)12 
(p 27)22 

(%f)h 

( Y 35v)^ 


1.17 x 10 16 
6.68 x 10 15 
-6.78 x 10 15 

4.12 x 10 11 
-1.84 x 10 16 

1.61 x 10 16 
5.27 x 10 15 
-7.07 x 10 13 

6.13 x 10 15 
-1.58 x 10~ 3 

-0.723 

0.703 

-0.676 

-0.371 

0.363 


-4.17 

X 

10 15 

3.98 

X 

10 15 

-4.90 

X 

10 15 

-5.52 

X 

10 12 

1.38 

X 

10 16 

1.49 

X 

10 16 

-1.69 

X 

10 16 

8.44 

X 

10 14 

-1.78 

X 

10 16 

1.50 

X 

10- 1 



1.93 


- 

1.19 


c 

1.730 


C 

1.733 


-C 

1.287 
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Table 3: Two example points in the parameter space presented in terms of the original 
parameters in the superpotential. 


parameter 

point 1 

point 2 

W& 351 ' [GeV] 

1.17 x 10 16 

-4.17 x 10 15 

m 2 7 [GeV] 

-2.92 x 10 14 

-1.64 x 10 15 

777.78 [GeV] 

-6.19 x 10 16 

-5.04 x 10 15 

Ai 

-6.49 x KT 1 

-8.91 x 10“ 2 

A'2 

5.66 x 10” 2 

-1.24 x 10” 1 

CO 

- 1.5 x nr 1 

2.78 x 10" 1 

A 4 

-7.38 x 10" 5 

-5.04 x 10" 4 

A 5 

-i.58 x nr 3 

1.50 x 10" 1 

^6 

-7.59 x KT 3 

3.06 x 10- 2 

A 7 

-4.23 x KT 1 

9.55 x 10" 1 

00 

5.08 

1.16 

(^ 27)11 

-7.23 x KT 1 

1.87 

(^ 27)12 

7.03 x KT 1 

-1.09 

(p 27)22 

-6.76 x 10" 1 

6.30 x 10- 1 


-3.71 x 10” 1 

7.39 x 10- 1 

&).. 

0 

0 

(^35l 7 )22 

3.63 x 10” 1 

-2.58 x 10" 1 
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Tabic 4: Tabic of predictions for two example point in parameter space. All charge fermion 
masses are in units of GeV, while the neutrino masses are in units eV. 


quantity 

experiment 

Vi 

parameter point 1 

fi(x) Xi 

parameter point 2 
fi(x) Xi 

rn c 

0.236 

0.226 

-0.432 

0.205 

-1.30 

m t 

92.2 

94.0 

+0.193 

105 

+1.38 

@cb 

0.0409 

0.0358 

-1.24 

0.0378 

-0.757 

m s 

0.013 

0.0144 

+1.06 

0.146 

+1.20 

rn h 

0.79 

0.79 

+0.0021 

0.781 

-0.112 

mv 

0.0599 

0.0613 

+0.241 

0.0664 

+1.08 

m T 

1.02 

0.867 

-1.51 

1.03 

+0.056 

m v3 


0.135 


0.0824 


m U2 


0.126 


0.0659 


(ml 3 - ml 2 ) 10 3 

2.32 

2.45 

+0.568 

2.24 

-0.325 

sin 2 023 

0.386 

0.343 

-1.12 

0.327 

-1.52 

x 2 



0.76 


1.03 


7 Discussion 


What we presented here is a model we believe is a good candidate for a minimal super¬ 
symmetric renormalizable Eq GUT. Let’s see why this model can be considered as more 
minimal than our previous candidate HU. Although the total number of degrees of freedom 
is now larger (78 > 27 + 27), what really counts are the number of multiplets and, even 
more important, the total number of free parameters. The number of parameters in the 
Higgs sector is now 11 complex minus 5 phases due to held redefinitions. Yukawa sector 
adds another 3 real and 6 complex parameters. Together with one real gauge coupling 
we have thus a total of 33 real parameters. This is 7 more than in the minimal SO(IO) 
[El mi ns], but still 16 less than even in the simplified truncated version of Eq in HU. 

Is there any possible low-energy signature of this model? It is often said that Eq could 
have possible light extra generations, coming from remnants of the three copies of the 
27p. Light vector-like fermions of the SM group could indeed emerge out of E 6 if they 
were associated with anomaly cancellation in a TeV-scale extra U(l). In such a setup, the 
extra U(l) would be a linear combination of the two U(l)s present in the rank 6 group Eg. 
Phenomenology of such TeV scale E 6 motivated models has been extensively studied, see 
for eg. |28j . In the present setup, however, there are no intermediate U(l)s and any light 
vector-like states would be accidental. In fact if we try to get such light states from our 
solution, we find them hard to obtain. The question is, does the matrix for (as an example) 


down quarks in eq. (78) allow 4 or more zero eigenvalues once we limit all u’s to zero? 


This can in principle be obtained either by putting C 2 , —> 0 (remember that our solution 
already has c\ = = 0), or by imposing a vanishing determinant constraint to Yukawa 
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matrices. The first case points towards an SU(5) invariant vacuum^ while the second one 
constrains the Yukawa parameters and so a worse fit to data is expected. So we conclude 
that such an extreme, albeit interesting situation is unlikely, at least in the given vacuum 
solution. Although we cannot make the same conclusions in general, any possibility of light 
states, if available, will occur due to fine-tuning in the superpotential parameters. 

Another possibility for having light states could be to have flat directions. We checked 
by explicit computation that no such states are present in our solution. Although we 
omit the details of this computation here, the interested reader can reconstruct the mass 
matrix of the SM singlets via the all-singlet terms in the superpotential, which are given 
in equations ([6])— (16). If the vacuum solution is then plugged-in, one discovers 4 massless 
singlet states, all of which are in fact would-be Goldstone bosons. Note that the adjoint 78 
of Eg contains 5 SM singlets, while only 1 remains unbroken among the SM generators, so 
the 4 massless singlet states are the would-be Goldstone bosons eaten up by the 4 broken 
singlet generators. We thus conclude there are no physical massless singlet states. Other 
vacua could in principle be possible: one of them is described in Appendix [Dj 

Proton decay is, as usual, quite hidden by details of superpartners’ spectra: it is hard 
to disentangle the GUT and SUSY breaking information from it. Once however colliders 
will (hopefully) tell us more about the low energy spectrum, this Eq theory as well as other 
grand unified theories could be tested better. For the sake of completeness the forms of 
the low-energy D = 5 operators are given in Appendix [E] 

The model we presented here is the minimal known E 6 , although for a more convincing 
proof we should satisfy three more checks. 


• First, there is a possibility to redefine the charges under matter parity, so that the 
27 + 27 parity is now transferred to the fermionic sector. This means that the fields 
351' + 351' + 78 should alone break to the SM gauge group as well as allow for DT 
splitting. The fermionic sector would now consists of 27 p, a = 1,... ,4 and 27 p. The 
Yukawa terms could be written schematically as 


27 f\ T ( y ^( 351 ') m 27 + A a (78)\ 

27 f) 1^27 + A 6 (78) i/asr'(351') ) 



(104) 


The total number of parameters is now the following: 5 complex parameters come 
from the Higgs superpotential, 3 phases of which can be rotated away by Higgs field 
redefinitions; 4 real diagonal components are given by the only Yukawa matrix, and 
8 complex and one real parameters are the off-diagonal terms; finally, 1 real gauge 
coupling sums to a total of 29 real free parameters. This would be 4 real parameters 
less than the model in this paper. Obviously there is no guarantee that such a model 
is realistic. We plan to come back to this issue in the future. 


• Second, we should study the complete three generation case, not only its two gener¬ 
ation subsystem. Although the number of parameters seems naively large enough, it 
is far from obvious that a successful fit is possible. In fact already in the two gen¬ 
eration case considered here we could not find a solution with vanishing y 1 2 , in spite 


1 This can be easily seen from the expression for A 2 in (971: if f 4 

both lead to an SU(5) invariant vacuum with e. 4 , —> 0, see (89). 


0 , we need either —> 0 or —> 0 ; 
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of enough free parameters. However the (at least partially) successful fit of the the 
minimal SO(IO) analogous case [29], [30] EH [32] [33] EH EH] EE] EE EE EH EE] EH] with 
10 and 126 (instead of 27 and 351') Yukawa couplings make us feel optimistic. Notice 
also that the usual obstruction of the neutrino mass either too low or unification 
violated ED- HU HE, IH] is here avoided as shown in section [6j 

• Third, the theory is not asymptotically free and has a huge gauge coupling beta 
function, more precisely 159. This means, similarly as in our previous Eq model, or 
the minimal renormalizable supersymmetric SO(IO), that a Landau pole is close to 
the GUT scale and so the theory itself may be already in the non-perturbative regime. 
An indication of problems being present already at the matching scale Mqut can be 
found in the calculation of the threshold corrections [ID] [HI] . The general problem of 
having large numbers of degrees of freedom (in our case through large representations) 
and the associated non-perturbativity is a problem, which is far from easy to solve 
and well beyond the purpose of this paper, although some progress has been made 
recently [12] BE] based on previous works on Seiberg dualities. We hope to come back 
to this very interesting issue soon. 

But even in the case the model presented here is the minimal one, other vacua could 
still be realistic with in principle different predictions. 
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A Gauge Boson masses 

We compute in this appendix the expression for the gauge boson masses, so that we can 
confirm that solutions really break into the SM group. We write the mass terms in the 
Lagrangian as 


C mass = g 2 A“M ab A» b , (105) 

where g is the Eg coupling constant and the matrix M ab is computed via 

M ah := ]TTr P<&») • (106) 

i 

The sum is over all representations containing VEVs (contributions come from the repre¬ 
sentations of the breaking sector), while t a denotes the action of the a-th generator on the 
representation 0j. The mass matrix becomes block diagonal if we choose an appropriate 
basis (indices a, b ), so that the basis gauge bosons have well defined transformation proper¬ 
ties under the SM group (note that some of these states are complex). We omit the details 
of this calculation and only give the results, collected in Table [5] 
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Tabic 5: Masses-squared of gauge bosons in SM representations using the ansatz 
ci = d,\ = e 2 = f 2 = e 5 = f 5 = ux = u 2 = y = 0. 


SO(IO) D SU(5) D 

SM D 

(mass) 2 /# 2 

45 

24 

(8,1, 0 ) 

0 

45 

24 

(1,3, 0) 

0 

45 

24 

(1,1, o) 

0 

45 

24 

(3,2,+f) 

KN 2 + i/ 4 | 2 ) 



(3,2,-f) 


45 

10 

(3,2,+i) 

>l 2 + !M 2 + !ld 2 | 2 


TO 

(3,2,-|) 

+ e ! 2 + A 2 + | e 4 2 + | / 4 2 

45 

10 

(3,1,-!) 

^H 2 + !lc 2 | 2 + !id 2 | 2 + | ei | 2 + i/ 1 | 2 


TO 

(3,1,+!) 


45 

10 

(1,1,+1) 

+H 2 + !| C2 | 2 + !ld 2 | 2 + | ei | 2 + l/ 1 | 2 


TO 

(1,1,-1) 


16 

10 

(3,2,+i) 

+ |\/I5w - v\ 2 + e 3 2 + |/ 3 | 2 + ! e 4 2 + 1|/ 4 | 2 

16 

TO 

(3, 2,-1) 


16 

10 

(3,1,-1) 

+ |\/15u; - v\ 2 + e 3 2 + |/ 3 | 2 + ! e 4 2 + ||/ 4 | 2 

16 

TO 

(3,1,+!) 


16 

10 

(1,1,+1) 

+ |\/l5u; - v\ 2 + e 3 2 + |/ 3 | 2 + 11 e 4 1 2 + ||/ 4 | 2 

16 

TO 

(1,1,-1) 


16 

5 

(3,1,+!) 

||u; + -^/3/5u 2 + ! \c 2 1 2 + | \d 2 2 + e 4 2 + \fi 2 + 

16 

5 

(3,1,-!) 

+ e 3~ + 1/3 2 + !M 2 +! 1/4 2 

16 

5 

(1,2,-!) 

il w +\/3/5 v 2 + !|c 2 | 2 + !|d 2 2 + |e 4 2 + |/i 2 + 

16 

5 

(1, 2, +!) 

+l e 3 1 2 +1 «/*3 1 2 + !M 2 + !i/ 4 | 2 

45 

1 

(1,1, 0) 

They mix: 

1 

1 

(1,1, 0) 

| [(A + B)±yj(A + By--%AB), 




A = 41 e 3 1 2 + 4 / 3 2 + e 4 2 + / 4 2 

B = 4 e 4 2 + 4 / 4 2 -|- c 2 2 + d 2 2 

16 

1 

(1,1, o) 

They mix: 

16 

1 

(1,1, o) 

!((C + D + \F\ 2 ) ± y/{C - Df + 16|i?| 2 ), 




C = |c 2 | 2 + 2 |e 4 1 2 + 2 / 3 | 2 + e 4 2 

D = \d 2 2 + 2|/i 2 + 2 e 3 2 + / 4 2 

E = eie 3 * + /T /3 
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B Particle identification 


Tabic 6: Identification of doublets and triplets in the representations of the Higgs sector. 


label 

E 6 D SO(IO) D SU(5) 

label 

Eg D SO(IO) 2 SU(5) 

doublet 

triplet 

Do,T 0 

78 D 16 D 5 

Dg, T 0 

78 D 16 D 5 

1 f 31 1 4 . 0 . 

+ 12 t “ ’ +12 31 

DuT, 

27 3 10 D 5 

Di,Ti 

27 3 10 D 5 

L' c 

d! 

D21T2 

27 D 10 D 5 

D2, Tg 

27 3 10 D 5 

V 

d ,c 

D 3 ,T 3 

27 D 16 D 5 

D3, T3 

27 3 163 5 

L 

d c 

L>4, 1 4 

351' 3 105 5 

D4, T 4 

35T 3 103 5 

< 0*0 

1 1 

e o H 

1 1 

a- H 

0 

Dg,T 5 

351' 3 103 5 

Dg, t 5 

351' 3 103 5 

Qu c — Lv c — AL's 
u c e c —d c v c +QQ—Ad' c s 

Dg,Tg 

351' 3 163 5 

Dg, Tg 

351' 3 163 5 

—Ls 

—d c s 

d 7 ,t 7 

351' 13 126 D 5 

d 7 , t 7 

351' D 126 D 5 

—Qu c - 3 Lu c 
—u c e c — 3 d c v c — QQ 

D 8 , T 8 

351' D 126 D 45 

d 8 ,t 8 

351' D 126 D 45 

Qd c + 3 Le c 

QL + u c d c 

Dg, Tg 

351' D 144 ^ 5 

Dg, Tg 

351' D 144 D 5 

—Qd' c + 4 L' c v c + L'e c 
—QL' + u c d' c + 4 d'v c 

fal 

O 

HI 

O 

351' 3 1443 5 

Dio, Tio 

351' D 144 D 5 

-L'v c 

-d! c v c 

-Du, Tn 

351' D 144 D 45 

Dn, Tn 

351' D 144 D 45 

—dd' c — 3e'e c 
-QL' - u c d' c 

T12 

351' D 126 D 50 

T\2 

351' D 126 D 50 

/ 

2u c e c - QQ 


In this appendix, we write the definitions of various states needed in the paper. The SM 
singlet VEVs were already defined in Table [l] We supplement the list of definitions with 
Tables [6] and [7j In Table [6] we define the doublets (1,2,+1/2), antidoublets (1,2,—1/2), 
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triplets (3,1, —1/3) and antitriplets (3,1, +1/3); these definitions are needed for DT split¬ 
ting. In Table [7j we define the weak triplets (1,3, ±1) relevant for type II seesaw. All the 
states in the two-index 351' can be specified by writing the basis states of this represen¬ 
tation by using two labels of the fundamental 27. In this notation, both labels are SM 
representation in the 27, assumed to have all the color and weak indices contracted in the 
correct manner to obtain the desired SM state in the 351'. More details on this notation 
can be found in [TTJ. 

Table 7: Induced VEVs in weak triplet scalars (1, 3, ±1) leading to seesaw type II. 


label 

E 6 D SO(IO) D SU(5) 

p.n. 

label 

E 6 D SO(IO) D SU(5) 

state 

Ax 

35T D 126 D 

15 

LL 

Ax 

351' D 126 D 15 

LL 

A2 

35T D 144 D 

15 

LL' 

A2 

351' D 144 D 15 

LU 

A3 

35T D 54 D 

15 

L'V 

A3 

351' D 54 D 15 

L'V 

a 4 

35T D 54 D 

15 

L ,C L ,C 

a 4 

351' D 54 D 15 

l' c l' c 


C DT splitting and Goldstone modes 

Analysis of the DT splitting in the group E 6 is complicated by the fact that a vacuum 
breaking to the SM group will automatically cause a doublet and a triplet mode to be 
massless. These massless modes are would-be Goldstone bosons: E 6 —>■ SM causes the 
breaking of 78 — 12 = 66 generators, with a doublet-antidoublet and triplet-antitriplet pair 
among them. One possible procedure to compute the condition for an extra massless mode 
in a matrix M, for which det M — 0, is to take 

lim £ ^o (det (M - el)/e) 

<e|/) ’ 1 J 

where e and / are the already present left and right null-eigenvectors of M, respectively. We 
present below, however, a simplified procedure of computing the conditions of DT splitting 
in the presence of a Goldstone mode. Its advantages are that it is computationally less 
intensive and that no issues with singularities, such as (e|/)=0, arise in the procedure. 

Suppose we use a generic label A for an n x n complex matrix. Although the true scalar 
mass-squared matrix is in fact the hermitian and positive definite matrix A^A, it is more 
efficient to work with A. A might not necessarily be diagonalizable, but it has a singular 
value decomposition. The presence of a zero eigenmode in A'A implies 

det A = 0. (108) 

We will rotate this matrix into a basis, where the left and right Goldstone modes corre¬ 
spond to the first basis vector of the rows and columns, respectively. We first write A in 
(1 + (n — 1)) x (1 + (n — 1)) block form: 

A=( m mR, V 

\m L M )' 
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(109) 










where boldface small letters denote n— 1 column vectors, and boldface capital letters denote 
matrices. Since A has a massless mode, there exist left and right null-eigenvectors e 0 and 
/o, respectively: 

A f e 0 = 0, Af 0 = 0. (110) 


Providing we choose the phases such that the first components of eb and /o are positive, 
we define (n — 1) column vectors e and f via 


£b 

e 0 




A = A/T^ftf 

l/ol A f 


(in) 


Writing the null-eigenvector conditions in equations (110) in block form, we get 3 indepen¬ 
dent equations (two (n — l)-vector, 1 scalar): 


A - ete mR + M' e = 0, (H2) 

\/l - ftf m L + M f = 0, (113) 

m\J\ — ete\/l — ftf - e'Mf = 0 (114) 


These can for example be used to define m, itll and mp in terms of M, e and f. We now 
define a n x n unitary matrix f/(x), whose form will be useful for rotating the basis of A: 
in block form, U is written as 


where 


tf(x) 


VT 


- x t x 

X 



A(x) 


I - 


xxi 

1 + a/ 1 — xfx 


(115) 


(116) 


It is possible to check explicitly that U is indeed unitary and that the following relations 
hold: 


f/(x)f/(x) t = J, 

(117) 

-1 = t/(x)t = U(-x), 

(118) 

A(x)t = A(x) = A(-x), 

(119) 

A(x) -1 — I + , _ xxt ,_ , 

y/1 — xtx (1 + \/1 — xtx) 

(120) 

detA(x) = 1- 

V ; l + Vl-xtx 

(121) 


The properly rotated matrix A, denoted by A ', is then by explicit computation equal to 


A! := U{e)AU{ f)- 1 


0 Ot \ 

0 M rot J ’ 


( 122 ) 


29 

















where 


M rot = A(c)-'M4(f)-‘ 


(123) 


We have indeed rotated into a basis, where the first column and row correspond to the 
zero eigenmodes. Crucially, the formula for the rotated (n — 1) x (n — 1) block has only 
one term, where a simple biunitary rotation is performed on the original block M; another 
zero-eigenmode can now be simply imposed by taking detM rot = 0. But since we now have 

detM = detM rot det(A(e)A(f)). (124) 


it is sufficient to impose det M = 0. The only possible caveat is the possibility that either 
detA(e) = 0 or detA(f) = 0; considering equation (121) and that e and f are parts of 
normalized vectors, this can happen only if ffie = 1 or f T f = 1, which would imply that the 
zero modes of A have a zero component in the direction of the first vector of the original 
basis. But since the eigenvalues of AA A (the determinant of A'A — XI) do not change if 
we rearrange the rows or the columns of A, we can always rearrange the original basis so 
that we take the z-th row and the j-th column to be the preferred one for the left and right 
null-eigenmodes, respectively. 

The main result thus states the following: if an n x n complex matrix A has a zero 
mode present, an additional zero mode is obtained by demanding the (i, j)-th minor of A to 
vanish (the subdeterminant of A, when the i-t\i row and j-th column are removed), where 
the indices i and j can be arbitrarily chosen, as long as the left and right null-eigenvectors 
of A have respectively a nonzero i-th and j-th component. 


D An alternative vacuum 


The Yukawa fit presented in this paper is based on the vacuum computed in section [3} Since 
we have not been able to obtain (yet) a full classification of all vacua in this E 6 model, 
we cannot conclude much about the broader possibilities of suitable vacua (especially the 
Yukawa sector), but we were able to End one other alternative Standard Model vacuum. 
We can obtain it by taking the following ansatz of vanishing VEVs: 

c 2 = d 2 = e 5 = f 5 = e 2 = h = «i = u 2 = y = 0. (125) 


This ansatz is similar to the ansatz of the original vacuum, but the vanishing of e^, fr, is now 
paired up with the vanishing of c 2 ,d 2 instead of c\,di. For the purposes of this appendix, 
we omit the specific form of the solutions, as well as other details, such as the check that 
the unbroken group is indeed that of the Standard Model. 

Looking at the Yukawa terms in equation (78), we see that the alternative vacuum with 
c 2 = h = 0 decouples the 16f from the 10f of SO(IO) (in leading order of m,Ew/M gut)- 
The heavy vector-like exotics (in the down-quark sector and charged lepton sector) in 
the fermionic 27f are thus purely in the 10f part. The low-energy mass matrices are thus 
simply those for the 16f, the analysis of the Yukawa sector becomes linear, and numerically 
one can make use of the fit for the minimal supersymmetric SO(IO) model [231, 30, |3I] 1321 
13311331135] 1361 H9| 120113711331139] with the Higgs in the 10 and 126 coupling to fermion pairs 
in 16 f- 

The intriguing possibility of this alternative vacuum, which recovers the SO(10) limit, 
is somewhat marred by issues in DT splitting. Given the ansatz in equation (125), the 
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doublet and triplet mass matrices become block diagonal with the following block form 
(with the basis of the barred states rearranged in the same order as for unbarred): 


Do Do Dr Do D 


T n T, Tr To T 


-Dio 

D 

Dio 

Dn 


n 


D i D 2 D 4 D* 


T To T a Tr 


d 7 D, 


t 7 n t 


12 


(126) 

(127) 


The Goldstone modes, for example, turn out to be in the first block. The problem now 
is, however, that a fine-tuning is block specific: the light Higgs lives only in one of the 
blocks. We see from equation (78) that the low-energy fermionic mass matrices M E and 
M e are now controlled solely by the EW VEVs v 2 , u 4 and b 8 . In the v s = 0 case, we get 
the unwanted mass relation = M E , but D 8 is in a separate block compared to D 2 and 
E 4 . A realistic pattern of fermion masses would thus require a double fine-tuning: one in 
the second block and one in the third block. Due to this feature, we consider this vacuum 
to be of less interest: beside the aesthetically unpleasing extra fine-tuning, the additional 
light Higgs pair H' u + H' d pair gives large threshold corrections in the running of the gauge 
couplings, possibly spoiling unification. The model contains many heavy states though, so 
the situation regarding the RGE is not clear-cut. For determining the feasibility of this 
vacuum, further investigation would be necessary. 

We conclude this section of the Appendix with a brief elaboration on which kind of 


ansatz is suitable for a good vacuum, i.e. we motivate equations (29)-(30) and (125). The 


considerations will be very similar to the ones in the E 6 breaking sector, where the 78 is 
omitted m and a full classification of vacua is known. Suppose we look for a specific 
vacuum solution: we want it to be as simple as possible (it has as many vanishing VEVs 
as possible), yet it needs to be able to break E 6 to the SM group. Due to the E-terms 
in SUSY, we assume a conjugate-symmetric ansatz, where a vanishing VEV in 27 or 351' 
implies that the corresponding (conjugate) VEV in the 27 or 351' also vanishes, and vice 
versa. First, we identify the SU(5) breaking VEVs from Table 0 e 4 , e 5 , / 4 , / 5 and y. The 
F y equation of motion automatically implies y = 0. Since SU(5) needs to be broken, either 
the pair e 4 , / 4 is non-vanishing, or the pair e 5 , f$. Considerations of alignment symmetry 
in m imply that the choice is irrelevant, since these pairs are exchanged if one changes 
the embedding of SO(IO) in E 6 , such that the two 5s of SU(5) in the 27 are exchanged. 
Therefore one pair needs to be necessarily non-zero, while we can try a vanishing ansatz 
for the other pair. In the 27 + 27 part of the breaking sector, the pairs c 4 , d\ and c 2 , d 2 
also get exchanged under alignment symmetry; we assume one pair to be non-zero (such 
that the 27 + 27 part does indeed contribute to the symmetry breaking), but we can again 
try setting the other pair to vanish in the simple ansatz. The choice of the vanishing pair 
now points to either the scenario of the main vacuum of this paper, or the the alternative 
vacuum presented in this appendix. The remaining part of the ansatz, e 2 — f 2 — 0 (which 
transforms into itself under alignment symmetry), is suggested from the solutions in the 
model without the 78, while u\ = u 2 = 0 then follows as a consequence of the E-terms. We 
finish the ansatz discussion with the following points: 


• The main vacuum solution of the paper follows from the ansatz, which is a direct 
extension of the solution ansatz in the absence of the 78 HU- There, the EOMs were 
simpler, and this ansatz actually represented the most general SM solution once gauge 
freedom and the E-term equations were accounted for. Note that only the ansatz can 
be extended, the main solution itself is not merely an extension of the solution when 
the 78 is omitted. 
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• The alternative ansatz leads to a SM vacuum only after the 78 was included; there 
is no such option if the 78 is omitted. 

• A complete classification would tell us, whether still other nonequivalent vacua exist, 
where the VEVs are non-vanishing. If the 78 is omitted, such vacua did not exist, but 
we expect this situation to change due to more terms and more VEVs in the EOM. 


E Proton decay 


For completeness let’s summarize the analysis of D = 5 proton decay in this model, similar 
to the analysis done in HU, obtaining analogous results. The low-energy operators in the 
superpotential, which are relevant for proton decay, are 


W | 


proton 


(CT A - C'; mA (X T ) m n ) [(1 + X*X T )-V\i (M^)ab 
(cf 4 + | c' 2 mjA (X T ) m n ) [(1 + IX*X T )- l / 2 ]J (Mr 1 ) 


C 


AB 


klB 


QiLjQ k Q i 


C ; 


klB 


d i u j u k e ii 

(128) 


with the A" defined in equation (85) and the mass matrix of M t as already defined in 
equation (55) of section [lj Note that the Mt has a zero eigenmode corresponding to the 
would-be Goldstone, so its inverse cannot be directly computed; we instead write M^ 1 = 
huiAf-s-oc^-MT + M e/ T ) -1 , with e and / being the left and right column-eigenvectors of 
Mt, respectively. 

The C coefficients are computed to be 
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2C'l jA 

2 C^ A 
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(130) 

(131) 

(132) 

(133) 

(134) 


We see that the C'-coefficients are the same as the coefficients in DU. 0 we cross out 
the contributions from the extra Yukawa term in that model. More specifically, notice that 
there are no A = 0 contributions from triplets/antitriplets in the new representation 78, 
since the 78 is not present in the Yukawa sector (it does not couple to two 27y’s). 

Once the model parameters are fit to the experimental values of the fermion masses 
and mixings, the four-fermion amplitude mediating proton decay is fixed and thus poten¬ 
tially dangerous, but it can always be suppressed by implementing a split supersymmetric 
scenario without changing any other conclusion. 
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